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Polarization and kinematics in Cygnus A 
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ABSTRACT 

From optical spectropolarimetry of Cygnus A we conclude that the scattering medium 
in the ionization cones in Cygnus A is moving outward at a speed of 170±34kms _1 , 
and that the required momentum can be supplied by the radiation pressure of an 
average quasar. Such a process could produce a structure resembling the observed 
ionization cones, which are thought to result from shadowing by a circumnuclear dust 
torus. We detect a polarized red wing in the [O III] emission lines arising from the 
central kiloparsec of Cygnus A. This wing is consistent with line emission created 
close to the boundary of the broad-line region. 

Key words: galaxies: active - galaxies individual: Cygnus A- infrared:galaxies - in- 
frared:ISM - quasars: general - ISM:dust 



1 INTRODUCTION 

Dust in active galaxies provides a screen that can both ob- 
scure the central regions from direct view, but also scat- 
ter nuclear radiation into our line of sight. By studying the 
scattering-induced polarization of the emission lines, some 
aspects o f the kinematics of th e scattering medium can be 
deduced ijWalsh fc Clegglll994t) . In this paper we study the 
polarization and kinematics of line emitting regions in the 
radio galaxy Cygnus A. 

Cygnus A is one of the brightest and nearest powerful 
double-lobed radio galaxies. Its proximity allows us to spa- 
tially resolve the two narrow-line regions on either side of 
the nucleus with sufficient signal-to-noise to study the po- 
larization of the emission lines. Hence, it is an ideal object 
to study the physics of the narrow-line region. From radio 
observations we know that the western side is oriented to- 
wards us, wi th the radio jet inclined at about 50° to our 
line of sight dCarilli fc BartheJll99rj) . We will refer to the 
western side as the near side, and to the eastern side as the 
far side. 

The op tical continuum p olarization has been studied 
in detail by lOgle et all Jl997l hereafter 097). They find a 
strong blue, polarized continuum in the west, consistent with 
the presence of a medium power central quasar. In the east 
the emission is less polarized, which they ascribe to dilu- 
tion by a young stellar population. Both in direct and in 
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polarized light they detect a broad component underlying 
the narrow Ha emission line. In the west they also detect 
polarized broad H/3, in agreement with the orientation of 
Cygnus A, and confirming the presence of a hidden quasar. 

Near- and mid-infrared images show b iconical struc- 
tures , resembling the optical ionization cones fl.Tackson et alJ 
1998). In the near-in frared, the edges of the cones are bright 
and highly polarized jTadhunter et alll999lh indicating that 
the walls of the ionization cones are dusty, and illuminated 
by th e central quasar. In a later paper (Tad hunter et all 
2000), it is shown that polarization is only seen along one 
side of the bicone. No explanation has yet been found for 
this asymmetry. Mid-infrared imaging at 10 and 18 /im con- 
firms the existence of a c onical shape, which is especially 
prominent on the far side l|Radomski et, alJ l2002h 

Ionization cones are common in active galaxies. The two 
physical processes that can generate low density cones are 
obscuration close to the central source as in young stellar ob- 
jects, and suspected to be present in most active galaxies, or 
a bi-polar outflow as detected in planetary nebulae. Obscu- 
ration by a dusty torus is postulated as p art of the unifica- 
tion s cheme for radio-loud active galaxies jUrrv fc Padovaml 
1995). However, outflow of matter has been detected in 
broad absorption line quasars, where gas moving away from 
the nucleus causes broad absorption lines in the optical and 
UV spectrum. Cygnus A show s outflowing gas at larg er dis- 
tance which emits [OIII] lines ( Tadhun ter et al.lll994l) . Out- 
flows can excavate funnels in the galaxy, in which the ISM 
will be ionized by the central source. They can be driven 
by either the radiation pressure of the central quasar, or a 
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nuclear starburst and do not necessarily require any obscu- 
ration to produce sharp shadows. Using existing deep optical 
spectropolarimetry of Cygnus A, we can determine aspects 
of the kinematics of the narrow- line region, and infer the 
process that dominates the creation of the ionization cones. 

For the remainder of this paper we will use 
Ho =65 km s _1 Mpc -1 and go =0.5, and a redshift for 
Cygnus A of z — 0.0561, determined from the redshift of 
the nuclear [OIII] doublet in our data. 



2 THE DATA REDUCTION 

For our analysis we will use the same data as used in the 097 
paper. They consist of deep spectropolarimetry of Cygnus A, 
obtained wit h the Low Resolution Imaging Spect rograph 
polarimeter jOke et al.lll995t iGoodrich et alJll995TI at the 
Keck telescope during the night of UT 1996 October 4-5. A 
longslit of 1" width was used, aligned with the radio axis at 
PA=101°. Two sets of data were obtained, one in good see- 
ing (0.6") of 4 x 900 seconds, and a second with worse seeing 
(1.0") of 4 x 1200 seconds. Both sets contain observations at 
four different angles of the half- wave plate (HWP), i.e. 0°, 
22.5°, 45° and 67.5°. The spectra cover 3500 to 8500 A in 
the observed frame, with a dispersion of 2.5 A per pixel and 
effective resolution of 10 A. 

Data reduction on the spectra was performed separately 
for each orientation of the HWP, splitting the usual ordi- 
nary and extra-ordinary rays. Standard IRAF routines were 
used for bias subtraction, flat field division, wavelength cal- 
ibration from NeAr arc spectra, and geometrical correction. 
The final 3cr error in both the wavelength calibration and 
the geometrical correction is 0.6 A, which corresponds to one 
fourth of the spectral pixel size. Three 0.85" apertures are 
traced in all separate 2-dimensional spectra, located 1.2 kpc 
east of the nucleus, at the nucleus itself, and at 1.2 kpc west 
of the nucleus. The width of the apertures corresponds to 
1 kpc in the host galaxy. The apertures used for sky subtrac- 
tion are 10" from the nuclear emission region on either side. 
The two sky apertures are averaged before subtraction. Af- 
ter the complete reduction the two data sets with different 
seeing are combined to obtain a better signal to noise in the 
emission lines. No galaxy subtraction is performed. 

After suitable binning, the 8 rays are combined to form 
the Stokes parameters /, Q and U, from which we compute 
the fractional polarization. We have co rrected for the Ricea n 
bias by using a Monte Carlo technique llFosburv et al.ll993l) . 
This method provides a distribution function of U, Q, P 
and 6p, and gives an estimate of th eir errors. The complete 
technique is described in detail bv IVernetl i200 if) . The re- 
sulting value of 8p is corrected for the instrumental offset 
and calibrated using observations of the polarized standard 
star HD 155528. 



3 RESULTS 

In this paper we will focus on the polarization of the emis- 
sion lines. For the results and analysis of the continuum po- 
larization and its different spectral components, we refer to 
the 097 paper. The strongest narrow lines were selected to 
calculate the fractional polarization, polarization angle and 
polarized flux over the profile. The general observed trend 
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Figure 1. Polarization in the [OIII] doublet. From top to bot- 
tom is shown the total flux in erg cm -2 s — 1 A - 1 , fractional po- 
larization, the polarization angle of the E-vector in the western 
aperture and the polarized flux in erg cm -2 s~ 1 A -1 . In the total 
flux spectrum the apertures are: west=solid line, nucleus= dotted 
line, east=dashed line. In the remaining three plots the apertures 
are coded as follows: triangles=east, squares=nucleus, filled cir- 
cles=west. The polarization angle is only plotted for the western 
aperture, the trend is similar in the other apertures. 

in all these lines is well illustrated by the [O III] doublet at 
4959,5007 A, shown in Fig. □ 

We observe significant polarization of P — 2 ± 1% in 
the emission lines in the east and west. The fractional po- 
larization decreases from the continuum value towards the 
line center. In the nucleus P = 1 ± 1%, consistent with no 
polarization. The value of P — 2± 1% is consistent with the 
value found by 097 of ~ 1.3%. In all apertures we observe a 
peak in polarization ~ 20 A redward of the line center. In the 
east and nucleus this causes a redshift of a few Angstroms 
in the polarized line flux. In the western aperture a rotation 
in polarization angle of ~ 20° is observed over the line. The 
same trend is seen in the other apertures, albeit with larger 
errors. 

The wavelength shift of the polarized line is particu- 
larly prominent in the nucleus, where the polarized line cor- 
responds to a red wing observed in the total flux line. This 
wing, as shown in Fig|21 is not seen in the other apertures, 
and remains undetected in other strong emission lines. How- 
ever, these are all significantly weaker, or blended in the red. 
Hence, the presence of a red wing as seen in the nuclear 
[O III1 doublet cannot be excluded. 
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Figure 2. Demonstration of the polarized wing on the nuclear 
[O III] 5007 A line. The top panel shows a best fit Gaussian with 
fixed central wavelength (line) to the observations (pluses), and 
in the middle panel the residual of the fit is shown. The bottom 
panel displays the polarized flux in the line. Units as in Fig. 



Deviations in the fractional polarization can be induced 
during the data-reduction process, especially by residual 
wavelength calibration errors between the o- and e-rays. To 
assess whether the observed red peak in the polarization is 
due to such an error, we constructed a simple model that re- 
produces the effects of line broadening and wavelength shifts 
in a single ray on the resulting polarization of a Gaussian 
emission line. The only way to reproduce a feature as ob- 
served in Cygnus A, is by changing the line width and off- 
setting the line center in two of the ordinary ray frames. 
This is impossible to reconcile with the fact that all data 
are reduced identically. In addition, the offset needs to be 
10 times larger than the achieved wavelength calibration er- 
ror. Therefore, we conclude that the observed signature is 
intrinsic. 



4 ANALYSIS OF THE LINE POLARIZATION 

As discussed above, the fractional polarization of the emis- 
sion lines shows a peak redward of the line center. Sim- 
ilar behavi our has been observe d in the planetary nebula 
NGC7027 iWalsh fc Clegdll994l) . where it is interpreted as 
an outward movement of the scattering medium. If the cen- 
tral quasar in Cygnus A can provide enough radiation pres- 
sure to push the matter in the ionization cones outward, this 
implies a similar situation can occur. 
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Treating scattering as instantaneous absorption and re- 
emission, it is easy to see that scattering of photons by a 
particle moving away radially from the source of emission 
will produce a Doppler shift. The observer sees the photon 
with a wavelength X obs , given by: 

Xobs =Ao[l+-(l-cos0)], 

c 

in which 6 is the angle between the direction of motion of 
the scattering particle and the line of sight to the observer, 
v is the particle velocity, Ao is the wavelength at which the 
photon is emitted, and assuming v <C c. 

The effect of the moving dust on the resulting polarized 
spectrum depends only on the angle of observation, 0. The 
situation is illustrated in Fig. [3] The term 1 + -(1 — cosf?) 
reaches a maximum of 1+ — for backward scattering. To the 
observer, this occurs only on the eastern side of the galaxy, 
so here a maximum shift is expected in the polarized line 
flux. For scattering through 90° the term 1 — cosf? equals 
1, which occurs on both sides. Therefore, it cannot generate 
an observed difference between the polarized emission lines 
from either side. In the case of forward scattering, the cosine 
term cancels, and the wavelength of the scattered emission 
is equal to the emitted wavelength. This is only possible on 
the western side. For scattering through angles close to 
or 180° there is little polarization. However, the scattered 
emission on the both sides arises from an integral over scat- 
tering angles between 90° and either 0° or 180°, which can 
result in the 3% polarization in the emission lines. 

By measuring the offset of the peak of the line in 
polarized flux on the east, the expansion velocity of the 
scattering medium can be obtained. From radio observations 
we kn ow that the angle of the jet is ~50° dCarilli fc Barthell 
1996), which corresponds to an angle of ~140° between 
the direction of motion of the scattering medium and 
our line of sight. The offset is measured to be 5±lA, 
which yields an outflow velocity of 170 ± 34 km s" 1 . This 
is very small compared to outflow velocities of several to 
tens of thousands kms -1 that are general ly measured in 
broad-absorption line quasars (|ElvisJ [2000) , and suggests 
that the scattering medium is unrelated to the high velocity 
clouds that produce the UV and X-ray absorption lines in 
such systems. 

To understand the nature of the driving force, we will 
estimate the expected velocity of dust clouds driven by the 
radiation pressure from the central source. A grain will be 
pushed outward if it experiences a radiation pressure larger 
than the force of gravity. The radiative acceleration on a 
single grain is given by: 

na 2 QL , . 

5rad = 5 (1) 

(see lLamers fc Cassinellil Il999t p. 165), in which Q is the 
radiation pressure efficiency, which we will assume to be 
unity, a the radius of the grain, L the total luminosity of 
the central source, r the distance, c the speed of light and 
irid the mass of a single grain. The grain mass is given by 
m d = |7ra 3 p s , and the specific density of the grain mate- 
rial p s = 3.6 gem -3 for silicate grains and 1.8 gem' 3 for 
carbonaceous grains. Classic interstellar dus t grains extend 
a ran ge of sizes from 0.001 /im to 10 fim (lLaor fc Draind 
1993). We assume an average grain radius of 0.1 /Ltm as in 
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Figure 3. Cartoon of the scattering geometry in Cygnus A. The 
solid arrows note the direction of the velocity of the dust clouds, 
the dashed arrows the direction of the light. East is to the left, 
west to the right, to observer is down. 



iRadomski et alJ <l2002l) . to balance the effects of larger ac- 
celeration for smaller grains, and lower abundance of the 
larger grains. We also assume pure silicate dust, adding 
carbonaceous grains to the dust will increase the radiative 
acceleration. Furthermore, we adopt a nuclear luminosity 
L = 10 45 4 ergs -1 ba sed on the recent X-ray observations 
bv lYoung et alJ i2002T) . At a distance of 1 kpc, the radiative 
acceleration in Cygnus A equals 5.1 xlO -6 cms -2 . The grav- 
ity acceleration exerted on the grain depends on the mass 
enclosed within the radius at which the grain sits: 



—j- = Mr (Mo) 



1.4 x 10" 17 cms 2 



assuming again a distance of 1 kpc, and M r being the mass 
within that radius in units of solar mass. The radiative ac- 
celeration exceeds the gravity if M r < 3 x 10 Mg. This 
is equivalent to the mass of a large spiral galaxy, and 
about two orders of magnitude larger than the mass of the 
super-massive black h ole in Cygnus A of ~ 2 x 10 9 Mg 
(Tadh unter et aljl2003T) . which is expected to dominate the 
central potential. Thus we can safely assume that the dust 
will be accelerated by radiation pressure. 

With M r <C 10 11 Mg we can express the acceleration of 



a dust cloud in terms of the central luminosity, dust-to-gas 
ratio (5) and distance from the central source: 



^cloud 



N d Tva 2 QL 



4 7rr 2 c 



C_ 



where Nd is the total number of dust grains in an optically 
thin cloud, given by 



N d = 
and 
C = 



35 Afcioud 
47ra 3 p s 



3SQL 
16ivap s c 



= 1.4 x 1CT 



in cgs units, for the values of the parameters adopted above. 
We use a Galactic dust-to-gas ratio of S = 0.01. 

The solution of the equation of motion gives the cloud 
velocity as a function of distance to the central source: 



v(r) 



1 - 



ro 



where ro is the distance from the nucleus where the dust 
cloud is formed, and the terminal velocity Voo is given by 
v/ 2C/ro. Adopting a cloud velocity of 170 km/s at a dis- 
tance of 1 kpc, we find that the cloud must have formed at 
ro = 0.8 kpc. This implies that the scattering medium orig- 
inates within the 1 kpc aperture. If it had been created in 
the nuclear regions and accelerated from there, its veloc- 
ity at 1 kpc would have been much higher than observed. 
Presumably the dust from the dense layer at the edge of 
the bicone is entrained in the turbulent flow in the bicone. 
When it enters the radiation field of the central source, it 
will experience the radiation pressure. The dust clouds will 
then move away from the nucleus, while scattering the light 
from the central source. This is consistent with the hypothe- 
sis derived from the infrared appearance of the bicones that 
their central regions have be en cleared from dust by radia- 
tion pressure, as proposed bv lTadhunter et al] J200(jt . 

Two basic assumptions are made. First of all, we as- 
sume the scattering cloud to be optically thin. If the cloud 
were optically thick, it would have to be created closer to 
the central source, since its effective scattering surface will 
be smaller, and the coupling of radiation to the dust parti- 
cles will be less efficient. The second basic assumption is that 
the drift velocity is small compared to the thermal speed of 
the gas, which implies that all momentum that is transfered 
from radiation to the dust cloud is shared by all dust and 
gas particles via internal collisions. This assumption is valid 
if the gas density is above a critical value that can be de- 
rived from eq. 7.44 in lLamers fc Cassinellil |l999) by setting 
the drift velocity equ al to the thermal speed. Adopting a 
temperature of 10 4 K <IOsterbrocklll978T) the thermal speed 
of the gas is of order lOkms -1 . This results in a lower limit 
for the electron density of the gas of 40 cm -3 . The typical 
observed ele ctron density in narrow-line reg ions lies around 
10 3-4 cm -3 iOsterbrockfll97gl : lFerlandlll981T) . several orders 
of magnitude above the critical lower limit, so this assump- 
tion is valid. 

The use of polarization shows unambiguously that an 
outflow must be present in the ionization cones. In case of an 
infall of matter, the absorbed emission will have a blueshift: 
^obs ~ Ao[l — -(1 + cosfl)]. On the side where 9 is small 
a blueshift of the line in polarized flux is expected, while 
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for large 9 no net effect should be visible. In the geometry 
of Cygnus A this would correspond to a blueshift in the 
eastern aperture. Using polarimetry of forbidden line profiles 
thus provides a unique method to measure infall and outflow 
velocities in any astronomical object. 



5 CONCLUSIONS 

Spectropolarimetry of (active) galaxies provides a powerful 
probe of the kinematics of the different components present. 
It can be used to uniquely determine outflow or infall of 
matter in any sufficiently polarized astronomical object. 

Large scale mass outflow occurs in Cygnus A, driven by 
the radiation pressure of the central source. The outflow ve- 
locity is measured to be 170±34kms~ 1 . This causes changes 
in the fractional polarization over the narrow emission lines, 
showing a shift of the line peak in polarized flux to the red. 
The polarized line flux from the east is redshifted, while on 
the west no shift is observed. This is only consistent with an 
outflow, not with infall of matter. 

A previously undetected redshifted nuclear component 
is found in [O III] , which is highly polarized. This compo- 
nent is consistent with line emission generated very close to 
the broad-line region, which is obscured from direct view. 
It is consistent with previous o bservations of anisotropic 
[O III] emission ijHes et alJ Il993h , and requires an obscur- 
ing medium of some sort in our line of sight to the inner 
narrow-line region. This nuclear component is not observed 
in any other narrow line, but its presence cannot be excluded 
due to blending or lower signal-to-noise in the other lines. 

The detection of an outward moving scattering medium 
in the ionization cones suggests that their shape could be de- 
termined by an outflowing wind driven by the radiation from 
the central quasar coupl ing to the dust. This is consi stent 
with previous results by iTadhunter et alJ <1999t l200oll . As 
a consequence, the detection of ionization cones in active 
galaxies cannot necessarily be used as evidence in favour of 
orientation-based unification, although it does not exclude 
the presence of an obscuring torus. The high obscuration 
of the central regions in Cygnus A and other radio galaxies 
could also be attributed to large scale dust. Detailed near- 
and mid-infrared observations, and modeling the emission 
of obscuring tori are required to give direct evidence of their 
presence. 
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